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Abstract 

Three-dimensional topological insulators arc characterized by insulating bulk state and metallic 
C surface state involving Dirac fermions that behave as massless relativistic particles. These Dirac 

fermions are responsible for achieving a number of novel and exotic quantum phenomena in the 

o 

CJ topological insulators and for their potential applications in spintronics and quantum computations. 

^ It is thus essential to understand the electron dynamics of the Dirac fermions, i.e., how they interact 

with other electrons, phonons and disorders. Here we report super-high resolution angle-resolved 
^ photoemission studies on the Dirac fermion dynamics in the prototypical Bi2(Te,Se)3 topological 

o 

^ insulators. We have directly revealed signatures of the electron-phonon coupling in these topological 

^ insulators and found that the electron-disorder interaction is the dominant factor in the scattering 

!• process. The Dirac fermion dynamics in Bi2(Te3_a;Sea;) topological insulators can be tuned by 

varying the composition, x, or by controlling the charge carriers. Our findings provide crucial 
information in understanding the electron dynamics of the Dirac fermions in topological insulators 
and in engineering their surface state for fundamental studies and potential applications. 
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Topological insulators represent a new state of matter that has an insulating bulk state 
but a conducting surface/edge state which hosts an unique spin texture and a linear disper- 
sion relation compatible with massless Dirac fermions[TH5]. In addition to their potential 
applications in spintronics and quantum computingjH [5], the topological insulators also 
provide a promising platform for realizing a number of novel and exotic quantum phe- 
nomena including quantum spin Hall effect [6], fractional charge and quantized currently, 
dynamical axion field[8], magnetic monopole[9], Majorana fermions[TOl [Tl] and quantum 
anomalous Hall effect [12]. Some of the phenomena have been observed experimentally in 
two-dimensional topological insulators, such as quantum spin Hall effect [13] and quantum 
Hall effect [HI in the HgTe quantum wells[6|. However, in the case of three-dimensional topo- 
logical insulators like Bi2(Te, Se)3 family [T5HT8] . in spite of intensive efforts from transport 
experiments [T91 - I23] . attempts to observe electronic transport signatures of the surface Dirac 
fermions have been hampered by a couple of obstacles. First, contrary to usual expectation 
that the bulk is insulating in topological insulators, the presence of defects in Bi2(Te, Se)3 
system makes the bulk conducting [23] and the transport properties are then dominated by 
the bulk state instead of the desired surface statefjOj. Second, in several cases where the 
surface transport signatures are identified, more than one frequency or one type of charge 
carriers contribute to the quantum oscillation [2T1423] . indicating the complexity of sample 
surface exposed to atmosphere or protection layers [25]. Third, the transport mobility of 
charge carriers in the Bi2(Te, Se)3 topological insulators, /i' ~ 10^ cm^V^^s^^'^^ I2S1 1^ . 
is not sufficiently high compared to those from HgTe quantum wells(1.5 x 10^ cm'^V^^s^^) 
[T3l [13] and graphene(2.3 x 10^ cm'^V^^s^^) [27I - I29] . Although strict electron backscat- 
tering is forbidden in topological insulators due to time reversal symmetry invar iance, it 
is important to uncover other scattering channels that limit the transport mobility of the 
two-dimensional Dirac fermions [26]. 

The charge transport properties are dictated by the electron djTiamics in materials. In 
particular, the scattering of electrons by other entities such as other electrons, phonons and 
disorder, determines the lifetime of the quasiparticles, and thus the electron mobility and 
mean free path in the transport process[30]. Angle-resolved photoemission spectroscopy 
(ARPES) has become a powerful tool to directly probe such many-body effects[3Tl 132] . 
Under the sudden approximation, ARPES measures a single particle spectral function 
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n [u — €k — ReT,{k,u)Y + [ImT,{k,u)]'^ 
with k and u representing electron momentum and energy, respectively. The real and imagi- 
nary parts of electron self-energy, ReTi{k, u) and Im'E{k, u), carry all the information about 
the band renormalization and quasiparticle lifetime in an interacting system pT| |22]- 

There have been a few investigations on the dynamics of Dirac fermions in Bi2(Te, Se)3 
topological insulators using ARPES[33l IM] or helium atom surface scattering probes [35]. 
However, the results obtained so far are not consistent. While the electron-phonon coupling 
strength (A) was reported to be 0.25 in Bi2Se3 from temperature-dependent photoemission 
measurements [33], it was claimed that the electron-phonon coupling is exceptionally weak 
in Bi2Se3 in another measurement [51] . The inconsistency comes presumably from two issues 
related to the ARPES study on many-body effects in topological insulators. First, the 
coexistence of bulk bands and surface bands near the Fermi level, due to the bulk doping 
by defects, complicates the identification of weak band renormalization effects in the surface 
bands. Second, the small momentum space and energy window of the Dirac fermions and 
the weak many-body effects require high precision ARPES measurements. The momentum 
space occupied by the surface state, ~ ± 0.1 A ^ for Bi2Se3 (Fig. lA) and ~ ± 0.04 A ^ for 
p-type Bi2Te3 (Fig. 2A), is rather small [TTl ES, [36] • To overcome these problems, it is ideal 
to first suppress the interference of the bulk bands in the measurements. Since the as-grown 
Bi2Se3 is electron-doped (n-type) while Bi2Te3 can be hole-doped (p-type) [25j , balancing the 
proportion of selenium and tellurium in Bi2(Se,Te)3 system may neutralize the bulk carrier 
doping and reduce the effect of the bulk bands on the surface state. Furthermore, vacuum 
ultra-violet (VUV) laser-based ARPES [37] can help suppress the bulk bands [25], in addition 
to its super-high energy and momentum resolutions (see Laser- ARPES measurements and 
Fig. S2 in Supplementary Materials) that are desirable in studying many-body effects in 
topological insulators. 

The Bi2Se3 sample shows sharp surface state band (Fig. lA) that facilitates the investi- 
gation of the Dirac fermion dynamics in the material. The momentum distribution curves 
(MDCs) at different binding energies (Fig. IB) show well-defined sharp peaks that can be 
fitted by Lorentzian lineshape to obtain peak position and peak width (FuU-Width-at-Half- 
Maximum, Tmdc) (see Fig. S3 in Supplementary Materials). The obtained MDC peak 
position as a function of the binding energy gives the measured dispersion relation (left 



panels of Fig. IC and ID). We assume a bare band as a linear line that connects the two 
points on the measured dispersion, one at the Fermi level and the other at 100 meV binding 
energy. This is reasonable within such a small energy window and particularly for the Dirac 
fermions that shows a linear dispersion over a large energy range. The effective real part 
of the electron self-energy (ReE) is obtained by subtracting the measured dispersion with 
the bare band, as shown in right panels of Fig. IC and Fig. ID (black empty circles). The 
imaginary part of the electron self-energy (ImS) is obtained from the fitted MDC width: 
ImE=rA^£)(7Xt>o/2 with vq being the Fermi velocity of the bare band (right panels of Fig. 
IC and Fig. ID, blue empty squares) [32] ■ 

The extraction of the electron self-energy provides full information for understanding the 
electron dynamics in a material where the electron scattering usually involves three main 
processes: electron-electron interaction, electron-phonon coupling, and electron-disorder 
interaction [381 ES]- The electron-disorder interaction is characterized by the imaginary 
part of the electron self-energy at zero energy (Fermi level) since the contribution from 
the electron-electron interaction and electron-phonon coupling approaches zero at the zero 
energy at very low temperature. The electron-phonon interaction is characterized by the 
band renormalization (usually a kink in dispersion, corresponding to a peak in the real part 
of the electron self-energy) and a drop in the imaginary part of the electron self-energy. 
As exemplified in Fig. 1 for Bi2Se3, there are clear indications of electron-phonon coupling 
revealed, such as the kink in dispersions (left panels of Fig. IC and Fig. ID), peak in 
the real part of the electron self-energy and a drop in the imaginary part of the electron 
self-energy (right panels of Fig. IC and Fig. ID). The characteristic energy scale of the 
phonons involved (wq) is determined by the peak position of the real part of the electron 
self-energy; in the case of Bi2Se3, it is ~18 meV as marked by arrows in Fig. IC and ID. 
The electron-phonon coupling strength. A, can be determined by the slope of the line near 
the zero energy in the real part of the electron self-energy. As shown by red dashed line in 
the right panel of Fig. IC, it is determined to be A~0.17 for Bi2Se3. The electron-electron 
interaction manifests itself more clearly in the imaginary part of the electron self-energy. In 
a typical Fermi liquid case, its contribution is proportional to the (E-Ei;')^: Imee=/3(E-Ei?)^. 
Since the electron-phonon contribution is constant above the characteristic phonon scale, the 
deviation from the constant above the phonon scale gives a measure of the electron-electron 
interaction involved. For Bi2Se3 and the Bi2(Se,Te)3 series, we find that the electron-electron 
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contribution is weak (/3~0.014 for Bi2Se3), particularly when compared with the dominant 
electron-disorder interaction and the electron-phonon coupling. Therefore, in the following, 
we will mainly focus on the electron-phonon coupling and electron-disorder interaction. 

As shown in Fig. 1, our measurements clearly indicate that there is an obvious electron- 
phonon coupling (A~0.17) in the electron dynamics of the Dirac fermions in Bi2Se3. It is 
nearly isotropic because similar behaviors are observed along the T-K (Fig. IC) and F-M 
(Fig. ID) momentum cuts. The characteristic energy scale, ~18 meV, is consistent with 
the highest-frequency phonon mode (^ig) at ~ 20 meV in Bi2Se3 from both theoretical 
calculations 110] and Raman measurements [IT] (see Fig. S7 in Supplementary Materials). 
This is in strong contrast to a recent report where exceptionally weak electron-phonon 
coupling was claimed and no dispersion kink was resolved in Bi2Se3[34J- We believe that 
the clear manifestation of electron-phonon coupling in our measurements is mainly due to 
much improved instrumental resolution we have achieved. 

The electron-phonon coupling is extremely weak in the p-type Bi2Te3 sample, as seen in 
Fig. 2. This sample was prepared under similar condition as the above Bi2Se3 sample and it 
is p-type because the Fermi level intersects with the bulk valence band[25]. Even measured 
with such a super-high instrumental resolution (see Fig. S2 in Supplementary Materials), 
we do not resolve any clear signature of electron-phonon coupling. The measured dispersion 
is basically a straight line, corresponding to a real part of electron self-energy that is nearly 
zero (black empty circles in right panels of Fig. 2C and Fig. 2D). There is no clear drop in 
the imaginary part of the electron self-energy (blue empty squares in right panels of Fig. 2C 
and Fig. 2D). All these observations point to an extremely weak electron-phonon coupling 
in the p-type Bi2Te3 sample. In the mean time, the p-type Bi2Te3 also exhibits a rather 
weak electron-disorder scattering when compared with Bi2Se3. As seen from the right panel 
of Fig. 2C, the imaginary part of the electron self-energy at zero energy is only 5.6 meV 
along the T-K direction. This is significantly lower than 37 meV along the same direction 
in Bi2Se3 topological insulator (right panel of Fig. IC). 

The dramatically different behaviors between Bi2Se3 (Fig. 1) and p-type Bi2Te3 (Fig. 
2) motivated us to further investigate the composition dependence of the Dirac fermion 
dynamics in the Bi2Te3_a;Sea; series. We observe a systematic evolution of band structure, 
electron-phonon coupling, and electron-disorder scattering, with the change of composition 
X in Bi2Te3_j;Se2: (Fig. 3). Note that all these samples are prepared under similar conditions. 
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First, the location of the Dirac point (E^) moves monotonically to higher binding energy 
with increasing x. For < a; < 0.05, the samples are p-type because the Fermi level intersects 
with the bulk valence band. For 0.15 < x < 1.5, the Fermi level intersects with the surface 
state bands only. For 1.5 < x < 3, the samples become n-type because the Fermi level 
intersects with the bulk conduction band. Note that the Dirac point exhibits a big jump 
when X changes from 0.05 to 0.15. Second, there is an overall increase of electron-disorder 
interaction with increasing x in Bi2Te3_2;Sei.. This can be seen from the increase of the 
imaginary part of the electron self-energy and hence its value at zero energy with increasing 
X. This indicates that the Bi2Te3_^Se2: system gets more disordered with increasing x. Third, 
there is also an overall increase in the electron-phonon coupling strength with increasing x 
in Bi2Te3_a;Sei:. Although no indication of electron-phonon coupling is observed in p-type 
Bi2Te3, a slight substitution of Te by Se, with x=0.05, already renders a clear observation of 
electron-phonon coupling in Bi2Te2.95Seo.05, as seen in the appearance of a peak in the real 
part of the electron self-energy (black circles in Fig. 3C2). The electron-phonon coupling 
strength increases from nearly zero for p-type Bi2Te3 to A~0.17 for Bi2Se3. 

It is interesting that for the same Bi2Te3 nominal composition, when prepared under 
different conditions, samples can become either p-type or n-type that exhibit dramatic dif- 
ference in their Dirac fermion dynamics. As shown in Fig. 3B, the Dirac point of the 
n-type Bi2Te3 lies nearly 310 meV below the Fermi level which differs significantly from 
the 85 meV location in the p-type Bi2Te3 (Fig. 3A1). Particularly, the n-type Bi2Te3 ex- 
hibits a pronounced electron-phonon coupling, as seen from the obvious peak in the real 
part of the electron self-energy (black circles in Fig. 3D) and a drop in the imaginary 
part of the electron self-energy (blue squares in Fig. 3D). The resultant electron-phonon 
coupling constant is A~0.19, which is even slightly larger than that in Bi2Se3 (A~0.17). 
The electron-disorder scattering also increases in the n-type Bi2Te3 compared to the p-type 
one, as the value of the imaginary part of the electron self energy at zero energy becomes 
much larger (Fig. 3D). However, compared with Bi2Se3, the disorder scattering in n-type 
Bi2Te3 remains much weaker(Fig. 3C11). The first principle calculation, which expects a 
rather weak electron-phonon coupling in Bi2Te3 (A~0.05)[42j, is hard to explain the strong 
electron-phonon coupling in n-type Bi2Te3 and particularly the disparity of electron-phonon 
coupling between the n-type and p-type Bi2Te3. As the main difference between the n-type 
and p-type Bi2Te3 lies in the relative position between the Dirac point and the bulk bands, 
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it is reasonable to deduce that the bulk bands may play a role in the electron dynamics of 
the Dirac fermions. This is consistent with the observation that, for n-type Bi2Te3 (x=3), al- 
though its composition is rather different from Bi2Teo.15Se2.85 (x=0.15), their Dirac fermion 
dynamics is quite similar in both the electron-phonon coupling and electron-disorder in- 
teraction (Figs. 3 and 4). The distinction of the electron-phonon coupling between the 
n-type and p-type Bi2Te3 indicates that the control of the charge carrier doping can be an 
alternative way to tune the Dirac fermion dynamics in topological insulators. 

Fig. 4 gives a quantitative summary of the electron-phonon coupling and electron-disorder 
scattering in the Bi2Te3_a;Sej; system. The relative position between the Dirac point and 
the bulk band is depicted in Fig. 4A for a couple of typical compositions and both p- and 
n-type Bi2Te3. With increasing x in Bi2Te3_a;Sea;, there is an abrupt jump in the electron- 
phonon coupling strength from x=0 to x=0.15, accompanied by a transition from p-type to 
n-type. Further increase of x leads to an overall slight increase of electron-phonon coupling 
(Fig. 4C). Meanwhile, the characteristic energy scale, Uq, also increases with x, which 
is consistent with the increase of the phonon frequency from Bi2Te3 to Bi2Se3 (see Fig. 
S7 in Supplementary Materials). This indicates that the (Te,Se) vibrations play a major 
role in dictating the electron-phonon coupling strength in the Bi2Te3_a;Sea; system. This 
is also consistent with the observation that the characteristic energy scale is identical in 
n-type and p-type Bi2Te3 because their phonon frequencies are quite similar (see Fig. S7 in 
Supplementary Materials) although they have different electron-phonon coupling strength. 
It is clear that electron-disorder scattering is dominant in the scattering process. As seen 
from Fig. 4B, the imaginary part of the electron self-energy at zero energy lmIl{Ep), which 
gives a good measure of the electron scattering rate, shows an overall increase from the 
p-type Bi2Te3 for x=0 to Bi2Se3 for x=3 in the Bi2Te3-xSex system (red circles in Fig. 
4B). The quasi-particle mean free path, which can be determined from the MDC width 
^MDC: — i/^MDC: fi^st dccrcascs abruptly from the p-type Bi2Te3 (x=0) to x=0.15, 
followed by an overall shght decrease with x in Bi2Te3_j;Te^ (blue circles in Fig. 3B). 
The electron mobility of the Dirac fermions (/x') is directly related to the quasi-particle 
mean free path, fi' = ei' /(hkp), with kp being the Fermi momentum. It is then estimated 
that the electron mobility in p-type Bi2Te3 is near 1000 cm'^V~^s~^ (see Supplementary 
Materials). By comparison, the electron mobility in Bi2Se3 and n-type Bi2Te3 is ~70 and 
~100 cm'^V~^s~^, respectively, which are significantly lower than that of p-type Bi2Te3. We 
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note that, although the electron mobihty obtained from ARPES is not strictly the same as 
the electron mobility in the transport properties, they have good correspondence, especially 
for materials with circular Fermi surface and nearly isotropic electron scattering. 

The present study will provide important information for engineering the Dirac fermion 
dynamics in Bi2Te3_a;Sea; topological insulators. It clearly indicates the presence of clcctron- 
phonon coupling and dominant role of the electron-disorder scattering in dictating the Dirac 
fermion dynamics in the system. The Dirac fermion dynamics can be tuned by varying 
either the composition or the charge carrier concentration. Specifically for the Bi2Te3_-rSea; 
system, our results clearly show that Bi2Te3 system is superior over Bi2Se3 in achieving 
higher electron mobility both due to its weak clcctron-phonon coupling and its weak electron- 
defect scattering. To further enhance the electron mobility in searching for novel quantum 
phenomena, it is necessary to choose a system with weak electron-phonon coupling and 
reduced disorder scattering. 
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FIG. 1: Observation of electron-phonon coupling in Bi2Se3. (A) Surface state dispersions along 
T-K (left panel) and T-M (right panel) directions. The inset shows the Fermi surface of Bi2Se3. 

(B) Momentum distribution curves (MDCs) of the surface state dispersion along T~K. Red line 
and blue line indicate the MDCs at the Dirac point (E/5) and the Fermi level (Ep), respectively. 

(C) Dispersion and effective electron self-energy of surface state band along T K direction. The 
left panel shows the measured dispersion (red circles) obtained from MDC fitting. The blue dotted 
line represents an empirical bare band that is a straight line connecting two points in the measured 
dispersion, one at Ej? and the other at 100 meV binding energy. The right panel shows the effective 
real part of the electron self energy, ReT,{k,uj) (black circles), and imaginary part of the electron 
self energy, ImT,{k,u}) (blue squares). The red dashed line is a linear fit to the low-energy region 
of ReT,{k,uj). Blue arrows indicate the characteristic energy scale of the electron-phonon coupling. 

(D) Dispersion and effective electron self-energy of surface state band along F-M direction. 
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FIG. 2: Extremely weak electron-phonon coupling in p-type Bi2Te3. (A) Surface state dispersions 
along T-K (left panel) and T~M (right panel) directions. The inset shows the Fermi surface of 
p-type Bi2Te3. (B) Momentum distribution curves (MDCs) of the surface state dispersion along 
T-K. Red line and blue line indicate the MDCs at the Dirac point (E^i) and the Fermi level (E^?), 
respectively. (C) Dispersion and effective electron self-energy of surface state band along T-K 
direction. The left panel shows the measured dispersion (red circles) obtained from MDC fitting. 
The blue dotted line represents an empirical bare band that is a straight line connecting two points 
in the measured dispersion, one at E^r and the other at 40 meV binding energy. The right panel 
shows the effective real part of the electron self energy, ReT,{k,u}) (black circles), and imaginary 
part of the electron self energy, Im'E,{k,uj) (blue squares). (D)Dispersion and effective electron 
self-energy of surface state band along F-M direction. 
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FIG. 3: Evolution of surface state band structure and electron self-energy with the composition, 
X, in Bi2Te3_a;Sex system. (A) The surface state band structure of Bi2Te3_a.Sei: with various 
compositions, x, measured along the T~K direction. (B) The surface state band structure of n- 
type Bi2Te3 along the T-K direction. (C) Corresponding effective real part of the electron self 
energy (black circles) and imaginary part of the electron self-energy (blue squares) for Bi2Te3_xSea; 
with different compositions. (D) Corresponding effective real part of the electron self energy (black 
circles) and imaginary part of the electron self-energy (blue squares) for n-type Bi2Te3. 
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FIG. 4: Variation of electron-phonon coupling and electron-disorder scattering with composition, 
X, in Bi2Te3_xSea; topological insulators. On the far right side, data on n-type Bi2Te3 are also in- 
cluded. (A) Schematic band structure at several typical compositions to show the relative position 
between the topological surface states (color dashed lines) and the bulk states (green solid lines). 
Black dashed lines indicate the Fermi level. (B) Variation of the electron-disorder scattering (red 
circles) with composition, x, in Bi2Te3_xSex topological insulators. It is obtained by taking the 
value of the imaginary part of the electron self energy at zero energy (Ef) as shown in Fig. 3C. 
Also shown is the electron mean free path (blue circles) calculated from the measured MDC width 
at the Fermi level. (C) Variation of the electron-phonon coupling strength A (red squares) and 
characteristic energy scale ujq with the composition, x, in Bi2Te3_2:Sea; topological insulators and 
in n-type Bi2Te3. The dashed lines in B and C are guide to the eyes. 
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